Unbounded growth of entanglement in models of many-body localization 
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An important and incompletely answered question is whether a closed quantum system of many 
interacting particles can be localized by disorder. The time evolution of an initially unentangled 
state is studied for a random-field XX Z Hamiltonian. Interactions induce a dramatic change in the 
propagation of entanglement and a smaller change in the propagation of particles. For even weak 
interactions, when the system is thought to be in a many-body localized phase, entanglement shows 
neither localized nor diffusive behavior but grows without limit in an infinite system: interactions act 
as a singular perturbation on the localized state with no interactions. The significance for proposed 
atomic experiments is that local measurements will show a large but non-thermal entropy in the 
many-body localized state, which develops over a diverging time scale as in glassy systems. 



One of the most remarkable predictions of quantum 
mechanics is that an arbitrarily weak random potential 
is sufficient to localize all energy eigenstates of a single 
particle moving in one dimension [HE]- In experiments 
on electronic systems, observation of localization is lim- 
ited to low temperatures because the interaction of an 
electron with its environment results in a loss of quan- 
tum coherence and a crossover to classical transport. Re- 
cent work has proposed that, if there are electron-electron 
interactions but the electronic system is isolated from 
other degrees of freedom (such as phonons), there can 
be a "many-body localization transition" even in a one- 
dimensional system for which all the single-particle states 
are localized [3HZ]- 

Progress in creating atomic systems where interactions 
between particles are strong but the overall many-atom 
system is highly phase coherent [5] suggests that this 
many-body localization transition may be observable in 
experiments [5]. Many-body localization is connected to 
the problem of thcrmalization in closed quantum systems 
as a system cannot thermalize if localization prevents 
equilibration of energy or particle number. 

The goal of the present work is to show that the many- 
body localized phase differs qualitatively, even for weak 
interactions, from the conventional, non-interacting lo- 
calized phase. The evolution of two quantities studied, 
the entanglement entropy and particle number fluctua- 
tions, show logarithmically slow evolution more charac- 
teristic of a glassy phase; however, the long-term behav- 
ior of these quantities is quite different. The growth 
of the entanglement entropy has previously been ob- 
served [101 111) to show roughly logarithmic evolution for 
smaller systems and stronger interactions. We seek here 
to study this behavior systematically over a wide range of 
time scales, showing that the logarithmic growth begins 
for arbitrarily weak interactions and does not saturate in 
the thermodynamic limit, and obtain additional quanti- 
ties that distinguish among possible mechanisms. We use 
matrix-product-state methods as well as exact diagonal- 
ization; the former allows some, not all, quantities to be 



computed for larger system sizes. 

One-dimensional s = | spin chains are a natural place 
to look for many-body localization [4 . To start, consider 
the XX model with random z-directed magnetic fields so 
that the total magnetization S z is conserved: 

ffo = J^£ (sfsf +1 + sfsf +1 ) + ]T hiSf. (i) 

i i 

Here the fields hi are drawn independently from the in- 
terval [—f],rj\. The eigenstates are equivalent via the 
Jordan- Wigner transformation to Slater determinants of 
free fermions with nearest-neighbor hopping and random 
on-site potentials; particle number in the fermionic rep- 
resentation is related to S z in the spin representation, 
so the z-directed magnetic field is essentially a random 
chemical potential. Now every single-fermion state is lo- 
calized by any 77 > 0, and the dynamics of this spin 
Hamiltonian are localized as well: a local disturbance 
at time t = propagates only to some finite distance 
(the localization length) as t —¥ 00. As an example, con- 
sider the evolution of a randomly chosen S z basis state. 
The coupling J± allows "particles" (up spins) to move, 
and entanglement entropy to develop, between two sub- 
regions A and B. But the total amount of entanglement 
entropy generated remains finite as t — > 00 (Fig. 1), and 
the fluctuations of particle number eventually saturate 
as well (sec below). The entanglement entropy for the 
pure state of the whole system is defined as the von Neu- 
mann entropy S = — tr p^ logp^ = — tr pg \ogps of the 
reduced density matrix of either subsystem. We always 
form the two bipartitions by dividing the system at the 
center bond. 

The type of evolution considered here can be viewed 
as a "global quench" in the language of Calabrese and 
Cardy 12] as the initial state is the ground state of an 
artificial Hamiltonian with local fields. Evolution from an 
initial product state with zero entanglement can be stud- 
ied efficiently via time-dependent matrix product state 
methods until a time where the entanglement becomes 
too large for a fixed matrix dimension. Since entangle- 
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ment cannot increase purely by local operations within 
each subsystem, its growth results only from propagation 
across the subsystem boundary even though there is no 
conserved current of entanglement. 

The hrst question we seek to answer is whether there is 
any qualitatively different behavior of physical quantities 
when a small interaction 



Hi n t — Jz SI S i+1 



(2) 



is added. With Heisenberg couplings between the spins 
(J z = J_l), the model is believed to have a dynami- 
cal transition as a function of the dimensionless disor- 
der strength f]/J z [4--6J. This transition is present in 
generic eigenstates of the system and hence exists at infi- 
nite temperature at some nonzero r\. The spin conductiv- 
ity or equivalently particle conductivity after the Jordan- 
Wigner transformation, is zero in the many-body local- 
ized phase and nonzero for small enough r)/J z . However, 
with exact diagonalization the system size is so limited 
that it has not been possible to estimate the location in 
the thermodynamic limit of the transition of eigenstates 
or conductivities. 

We find that entanglement growth shows a qualitative 
change in behavior at infinitesimal J z . Instead of the ex- 
pected behavior that a small interaction strength leads 
to a small delay in saturation and a small increase in 
final entanglement, we find that the increase of entan- 
glement continues to times orders of magnitude larger 
than the initial localization time in the J z = case (Fig. 
1). This slow growth of entanglement is consistent with 
prior observations for shorter times and larger interac- 
tions J z — 0.5Jl and J z = Jj_ [TU1 HT] , although the 
saturation behavior was unclear. Note that observing 
a sudden effect of turning on interactions requires large 
systems, as a small change in the Hamiltonian applied 
to the same initial state will take a long time to affect 
the behavior significantly. We next explain briefly the 
methods enabling large systems to be studied. 

To simulate the quench, we use the time evolving block 
decimation (TEBD) [T3 [H] method which provides an 
efficient method to perform a time evolution of quantum 
states, \ip(t)) — U(t)\tp(0)), in one-dimensional systems. 
The TEBD algorithm can be seen as a descendant of 
the density matrix renormalization group [15] method 
and is based on a matrix product state (MPS) repre- 
sentation [TBI [17] of the wavefunctions. We use a second- 
order Trotter decomposition of the short time propagator 
U (At) = cxp(-iAtH) into a product of term which acts 
only on two nearest-neighbor sites (two-site gates) . After 
each application, the dimension of the MPS increases. To 
avoid an uncontrolled growth of the matrix-dimensions, 
the MPS is truncated by keeping only the states which 
have the largest weight in a Schmidt decomposition. 

In order to control the error, we check that the ne- 
glected weight after each step is small (< 1CP 6 ). Al- 




FIG. 1. (a) Entanglement growth after a quench starting 
from a site factorized S z eigenstate for different interaction 
strengths J z (we consider a bipartition into two half chains of 
equal size). All data is for 77 = 5 and L = 10, except for J z = 
0.1 where L = 20 is shown for comparison. The inset shows 
the same data but with a rescaled time axis and subtracted 
J z — values, (b) Saturation values of the entanglement 
entropy as a function of L for different interaction strengths 
Jz. 



gorithms of this type are efficient because they exploit 
the fact that the ground-state wave functions are only 
slightly entangled which allows for an efficient truncation. 
Generally the entanglement grows linearly as a function 
of time which allows only the simulation of rather short 
times. In our case, however, the entanglement growth 
is logarithmic and thus it is possible to perform a time 
evolution over long time scales [11) . 

In order to study the asymptotic (t — > 00) behavior, we 
employ exact diagonalization techniques. We diagonalize 
the full Hamiltonian and construct U (t) which allows us 
to access very late times after the quench. Since the 
Hilbert space grows exponentially with the system size, 
we can only consider rather small systems. All the data in 
this paper are obtained by averaging over more than 10 4 
field configurations {hi} starting from a random product 
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FIG. 2. (a) Growth of the particle number fluctuations of 
a half chain after a quench starting from a site factorized 
5 Z eigenstate for different interaction strength J z . All data 
is for r\ — 5 and L = 10 except for J 2 = 0.1 where L = 
20 is shown for comparison. The inset shows the same data 
but with a rescaled time axis and subtracted J z — values, 
(b) Saturation values of the particle number fluctuations as a 
function of L for different interaction strengths J z . 



state 1-0(0)) = |mi)|rn 2 ) . . . \m L ), where \m,j) £ {| f),| !>} 
are eigenstates of the S*| operators (for the saturation 
data we always start from a Neel state to reduce noise 
- the data is qualitatively the same for a random initial 
state). 

In Fig. [T] upper panel, we plot the time evolution of 
the half-chain entanglement entropy. The initial quick 
rise of the entanglement till time J±t ~ 1 corresponds 
to expansion of wavepackets to a size of the order of 
the localization length. This rise is independent of the 
interaction strength J z . After this initial expansion, 
the entanglement saturates in the non-interacting case 
but increases logarithmically for any non-zero interac- 
tion strength. The time at which the logarithmic growth 
starts is J z t ~ 1, and all curves for different interacting 
strengths roughly collapse on a single curve when plot- 
ted against J z t (see inset). In the case of J z /J± = 0.1 



we have plotted the entanglement for both L = 10 and 
L = 20. While absolute value of the entanglement dif- 
fers slightly, the two curves are indistinguishable after 
subtracting the non- interacting value (inset), suggesting 
that for the entanglement at the time scales explored, 
finite size effects are small. 

In an infinite system, we expect the logarithmic growth 
of the entanglement to continue indefinitely. In a finite 
system, in contrast, it will saturate at late times at some 
value Soo. The lower panel in Fig.[T]shows the system size 
dependence of this saturation value, as obtained by exact 
diagonalization. In the non-interacting case, this value is 
independent of system size as expected if the localization 
length ^«L. In the presence of interaction the satura- 
tion value is independent of the interaction strength (we 
have checked that the small difference between the curves 
at large L goes away as the saturation time is increased) 
and follows a volume law. For small sizes the result is 
roughly consistent with Ref. [101 although going to longer 
times we find that the saturation follows a volume law 
rather than being logarithmic in subsystem size as con- 
jectured there. Such logarithmic dependence occurs at 
certain random-singlet quantum critical points in one di- 
mension [T5], but the present system is not critical by 
standard definitions. 

The observation of a volume law for the saturation 
entanglement entropy suggests at first glance that the 
system may be partially thermalized, as a finite entan- 
glement entropy per site is consistent with a thermal en- 
tropy. However, the saturation value of the entangle- 
ment entropy per site is much lower than the infinite- 
temperature thermal value associated with 2 states per 
site. Since, as shown below, particle number propagates 
much less than entanglement and may well be localized, 
the entanglement entropy at saturation may result from 
an incomplete thermalization where the system locally 
equilibrates among a restricted ensemble with particular 
values of particle number, energy, and other conserved 
quantities. This would be a many-body-localized version 
of the generalized Gibbs ensemble describing thermaliza- 
tion in some systems without disorder. 

The time evolution of the half-chain particle number 
fluctuations, i.e., the variance of the total spin on half the 
chain, is plotted in the upper panel of Fig. [2j The be- 
havior of the particle number fluctuations is qualitatively 
different than the entanglement entropy. Interactions do 
enhance the particle number fluctuations, but while there 
are signs of a logarithmic growth as for the entanglement, 
this growth slows down with time. The saturation values, 
similarly, are interaction dependent and only fractionally 
larger than in the non-interacting case. This does not 
seem to be a case of the system not having reached satu- 
ration, as increasing the saturation time by two orders of 
magnitude has a negligible effect on the saturation val- 
ues. The half-chain energy variance behaves in a man- 
ner similar to the particle variance (data not shown). It 
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would therefore seem that particle transport is not fully 
ergodic, which would be consistent with a many-body 
localized phase. 

Our results are consistent with a phase diagram con- 
sisting of a conventionally localized phase at zero inter- 
action, a quasi-localized or many-body localized phase 
for moderate interactions, and a delocalized phase when 
interactions dominate disorder. Since all local measure- 
ments on a subsystem are entirely determined by the re- 
duced density matrix of that subsystem, the entangle- 
ment entropy difference between Anderson and many- 
body localized phases will appear in any local measure- 
ment of entropy. Simple models can be constructed of a 
many-body localized phase with zero particle number and 
energy conductivities, but with infinite saturation entan- 
glement, compared to finite saturation entanglement in 
the non-interacting case |19j . 

We comment briefly on the possibility that the parti- 
cle and energy variances are not actually saturating but 
increase without limit, which cannot strictly be ruled out 
from our data. (This would still be consistent with zero 
conductivity if, for example, the continuing increase re- 
sults from subdiffusive motion of particles.) It should be 
noted that the quantities studied here by matrix-product- 
state methods are not directly comparable to quantities 
previously studied by exact diagonalization, e.g., the level 
statistics of eigenstates [3J [TU] or correlations in a single 
eigenstate [5]. It is conceivable that there is a transition 
as a function of rj in such quantities but that the real-time 
evolution studied here is insensitive to that transition. 
However, if this is the case, it means that experiments 
must be carefully designed to detect the transition as the 
dynamics of a simple initial state are characterized better 
by glassy dynamics than ordinary localization. 

If the ultimate dynamics on long time scales is de- 
phased, as one part of the system interacts with a large 
region that can serve as a "bath" , then classical dynam- 
ics is appropriate, and in classical physics it has been 
known for many years that quenched random systems 
need not be diffusive, especially in ID. One model lead- 
ing to anomalous long-time dynamics was introduced by 
Sinai [50], who showed that one-dimensional Brownian 
motion is modified by the addition of a static Gaussian- 
correlated random force F, 

d ^ = v (t) + F{x). (3) 

The root-mean-squared displacement R — \J (x 2 ) of a 
particle from its initial location scales asymptotically as 
R ~ log 2 (£) rather than R ~ \ft as in ordinary diffusion. 
The scaling of R can be obtained from a simple estimate 
of the time required to cross an occasional region of un- 
favorable potential [21 . Another possibly relevant clas- 
sical model follows from considering classical rate equa- 
tions with a broad distribution of rates [55] , which leads 
to logarithmic diffusion R ~ log t for intermediate times 



before ultimately crossing over to diffusion. 

In conclusion, the many-body localized phase is fun- 
damentally different in at least one measurable dynam- 
ical property, the creation of entanglement, from non- 
interacting Anderson localization. Two important ques- 
tions for future work are a full characterization of other 
dynamical properties of this phase, including whether 
particles and energy are fully localized or propagate sub- 
diffusively, and an understanding of how these proper- 
ties connect to the unusual dynamical phase transition 
between the many-body localized phase and the delocal- 
ized phase. The authors acknowledge conversations with 
P. Calabrese, D. Huse, V. Oganesyan, and G. Refael, and 
support from DOE BES DMSE (J. H. B) and ARO via 
the DARPA OLE program (J. E. M.). 
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